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Introduction 


In  the  latter  half  of  the  1960 's  the  field  of  manpower  planning 
experienced  a marked  shift  with  the  development  of  goal  programming 
models.-^  This  shift  is  documented  by  two  keystone  papers,  i.e.,  "A 
Goal  Programming  Model  for  Manpower  Planning"  by  A.  Charnes,  W.W.  Cooper  and 
R.  J.  Niehaus  /18.27,  and  "Static  and  Dynamic  Assignment  Models  with  Multiple 
Objectives  and  Some  Remarks  on  Organization  Design"  by  A.  Charnes, 

W.  W.  Cooper,  R.  J.  Niehaus  and  A.  Stedry  /l8.£7.  Since  the  1960's,  a 
continuing  series  of  papers  has  documented  the  extensions  and  follow-on 
efforts  which  are  still  in  progress  today!  The  purpose  of  this  paper  is 
to  review  these  developments,  discuss  the  current  status  and  project  the 
research  frontiers  which  appear  to  be  on  the  horizon. 

Before  moving  into  the  subject  of  this  paper,  it  is  necessary  to 
comment  on  the  people  by  whom  the  research  was  actually  accomplished.  In 
a very  true  sense  this  paper  could  be  authored  by  A.  Charnes,  W.  W.  Cooper 
and  R.  J.  Niehaus  rather  than  just  the  present  author.  The  work  has  been 
a team  effort  from  the  beginning.  A.  Charnes,  W.  W.  Cooper  and  R.  J. 

Niehaus  have  been  the  core  of  the  team,  which  has  been  augmented  and 
extended  by  others  who  participated  on  a full  contributing  basis  as 
particular  areas  of  emphasis  were  encountered.  Among  these  are  D.  Sholtz, 

G.  L.  Thompson,  A.  Stedry,  B.  Moore,  D.  Klingman,  L.  Mannis,  E.  Bres, 

D.  Cass,  K.  A.  Lewis,  A.  Nelson,  A.  Albanese,  K.  Padalino,  S.  Korn, 

E.  Leader,  B.  Hall  and  D.  Nitterhouse.  As  long  as  there  is  a clear 

1 J See  Charnes  and  Cooper  £107  for  a recent  summary  of  field  of  goal 
programming. 


understanding  that  the  work  reported  is  that  of  a large  complex  team, 
the  advantage  of  being  the  sole  author  is  that  the  work  can  be  discussed 
from  the  vantage  point  of  an  insider  and  at  the  same  time  express  my 
privilege  of  being  a colleague  of  Abe  Charnes  and  Bill  Cooper. 

There  are  two  types  of  developments  which  have  resulted  from  the 
Navy  Civilian  Manpower  Modeling  Program.  The  first  is  the  models  which 
have  been  subjected  to  extensive  prototypes  and  eventually  selected  for 
implementation.  The  second  is  the  models  and  ideas  which  have  been 
reported  in  the  literature  with  little  further  attention  paid  to  them. 

The  first  set  of  models  will  be  reviewed  briefly  with  the  appropriate 
references  so  that  thrust  of  the  implementation  is  clear.  Attention  will 
also  be  paid  to  the  second  set  since  many  good  ideas  have  been  proposed 
which  may  very  well  provide  the  linking  pins  to  future  developments. 

These  parallel  tracks  will  then  be  used  to  highlight  some  of  the  possible 
research  frontiers. 

II . Aggregate  Model  Structures 

The  initial  model  developed  by  Charnes,  Cooper  and  Niehaus  ^18. £7  was 
a goal  programming  model  with  embedded  Markov  processes.  The  variables 
are  defined  as  : 

E^(t) , E£(t)  * positive  or 'negative  deviation,  respectively, 
for  kch  manpower  category  on  time  t. 

Pkt  - the  weight  which  is  applicable  to  kch  manpower  category 
in  period  t. 

ni(t)  - the  ith  component  of  n(t)  where  n(t)  = accumulated 
inventory  from  hiring  to  period  t. 


g (t)  = net  requirement  for  the  kC  manpower  category  after 

allowance  for  the  available  persons  remaining  in  this 
category  growth  and  attrition  which  occurred  during 
the  t previous  periods. 

c^(t)  = the  ic^  mean  salary  value  in  period  t. 

B(t)  = applicable  budgetary  ceiling  for  salary  expenditures  in 
period  t. 

Mc  * Markov  transition  matrix  depicting  the  transfer  probabilities 
between  job  categories  where  Mt  = MM. . .M  (t  times) 

a = vector  of  initial  inventory  of  personnel  in  all  job 
categories. 

This  initial  model  can  be  stated  in  its  transformed  and  reduced 
version  as: 

(Objective  - Minimize  weighted  sum  of  deviations  from  net  manpower  goals:) 

Min  £ Z pkt  $+(t)  + E£(t)J 

Subject  to: 

(Goal  Constraints) 

Z nt(t)  - E+(t)  + E”(t)  = g,  (t) 
ielk  1 k k k 

(Manpower  Transition  Conditions) 

-(M)^  n(t)  + n1(t+l)  ^ 0 
(Salary  Budget  Constraints) 

Z Ci(t)  n±(t)  ^ B(t)  - cT(t)  Mta 
(Non-Negativity  Constraints) 

n^t),  E^(t),  E£(t)  ^ 0 

The  above  model  was  tested  via  numerical  examples  $37,  $8.37  and  a 
number  of  extensions  appeared  warrented.  The  first  of  these  involved 
revising  the  transition  matrix  to  Include  an  adjustment  for  retirements  $8.37 
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This  was  necessary  since  the  Navy  civilian  workforce  at  that  time  had  a 


large  age  bulge  of  personnel  nearing  retirement.  An  extension  was  also 
made  by  Charnes , Cooper  and  Niehaus  £8.47  to  include  training  explicitlv 
in  the  models.  Subsequently,  upon  closer  examination,  it  was  determined 
that  the  data  was  not  available  in  any  form  which  would  make  this  extension 
useful  to  studying  training  problems. 

The  training  extension  turned  out  to  be  highly  important  for  another 
reason.  The  initial  model  requires- a large  number  of  pre-computed  numbers 
(i.e.,  net  manpower  requirements,  adjusted  budgetary  numbers  etc.).  This 
means  that  time-consuming  and  thus  expensive  computer  processing  was 
necessary  for  pre-processing.  Even  more  relevant  was  the  fact  that  all 
these  intermediate  numbers  required  by  the  model  made  it  necessary  to 
provide  extensive  explanations  to  Navy  managers  with  a resulting  high 
level  of  confusion.  The  revised  generalized  network  model  permits  relevant 
numbers  in  terms  that  managers  can  understand  to  be  input  directly  without 
change.  Similarly,  there  was  no  need  for  translation  of  the  outputs.  In 
this  form,  the  net  numbers  previously  pre-computed  are  automatically 
developed  in  the  solution  of  the  linear  program. 

Another  missing  element  in  the  initial  model  was  the  ability  to  ac- 
commodate fires  or  Reductions-In-Force  (RIF's)  in  the  model.  This  became 
particularly  important  as  the  Navy  drewdown  its  civilian  workforce  in  the 
early  1970's  (from  425,000  to  300,000  in  three  years).  Also,  in  the  use 
of  the  model  for  skills  balancing,  particularly  in  the  naval  industrial 
facilities,  information  on  such  possible  adverse  impacts  is  desirable. 

The  generalized  network  form  lent  itself  to  this  change  since  hiring  and 
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firing  appear  explicitly  as  variables  which  can  be  appropriately  weighted 

in  the  objective  function.  This  version  can  be  found  in  Charnes,  Cooper, 

Niehaus  and  Sholtz  /T.8.57  and  in  Niehaus  £[8.l7.  The  Recruiting  Requirement 

Model  (RRM) , as  this  model  was  named,  has  become  the  workhorse  of  the 

operational  modeling  system.  The  variables  are  defined  as: 

Ej£(t),  E£(t)  = positive  or  negative  deviation,  respectively, 
for  kc^  manpower  category  in  time  t. 

x^(t)  = manpower  on-board  (in  place)  in  the  kth  manpower 
category  in  period  t. 

yjt(t)  = hires  in  kfc^  manpower  category  in  period  t. 

zfc(t)  = fires  or  Reductions-In-Force  (RIF's)  in  the  k^ 
manpower  category  in  period  t. 

Gjt(t)  = manpower  requirement  (goal)  for  the  kth  manpower 
category  in  period  t. 

akt  = t*ie  weight  which  is  applicable  to  the  positive  goal 

discrepancy  for  the  kth  manpower  category  in  period  t. 

= the  weight  which  is  applicable  to  the  negative  goal  dis- 
crepancy for  the  k*1^  manpower  category  in  period  t. 

^kt  = t*le  wei8ht  on  the  hires  in  the  k*-*1  manpower  category 
in  period  t. 

6^  = the  weight  on  the  first  or  Reductions  in  Force  (RIF's) 
in  the  kc^  manpower  category  in  period  t. 

M = Markov-like  transition  matrix  depicting  the  probabilities 
of  internal  personnel  movement 

a = vector  of  initial  inventory  of  personnel  in  all  job 
categories 

T T 

e = sum  vector  (i.e.  e = (1,1, ...1) 

C(t)  = manpower  ceiling  for  period  t. 

T 

s = average  salary  vector  in  period  t. 

B(t)  = salary  budgetary  ceiling  for  period  t. 
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The  structure  of  the  model  in  its  transformed  and  reduced  form  is 
as  follows : 

(Objective  - Minimize  weighted  sum  of  deviations  from  gross  manpower 
goals  and  weighted  number  of  hires  and  fires  (Reductions- 
In-Force)) 

Min  E E /akt  E+(t)  + 6kt  Ek(t)  + ykt  yk(t)  + 6kt 
Subject  to: 

(Goal  Constraints; 

*k(t)  - E+(t)  + Ek(t)  = (^(t) 

(Manpower  Transition  Conditions) 

xk(0)  = a 

- M xk(t-l)  + X^t)  - yk(t)  + zk(t)  = 0 

(Manpower  Ceiling  Constraints) 

eT  xk(t)  ^ C(t) 

(Salary  Budget  Constraints) 

sT(t)  ^(t)  ± B(t) 

(Non-negativity  Constraints) 

xk(t)’  yk(t)«  zk(t>’  Ek(t)»  Ek(t)  " 0 

This  basic  model  was  also  extended  by  Charnes,  Cooper,  Niehaus  and 

Sholtz  in  a number  of  ways.  One  was  the  extension  to  incorporate 

input-output  considerations  to  study  program  planning.  This  extension  was 
discussed  in  detail  by  Niehaus  /l8.£J , including  the  development  of  a 
numerical  prototype.  Work  progressed  through  the  programming  of  the 
necessary  extensions  to  the  linear  programming  matrix  generator  and 
report  writer.  Implementation  was  not  carried  beyond  this  point,  however, 
due  to  lack  of  resources  to  continue  its  development.  The  additional 
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variables  of  this  revised  model  are: 

(t) , (t)  = positive  or  negative  deviation,  respectively, 

for  the  i*^  final  user  support  category  in 
period  t. 


Vi(t) 

= 

amount  of  support  for  the  ith  final  user  in  terms 
of  its  dollar  budget  in  period  t. 

Wj  (t) 

= 

amount  of  output  of  the  jth  producer  in  terns 
of  its  dollar  budget  in  period  t. 

*i 

= 

the  weight  which  is  applicable  to  the  positive  goal 
discrepancy  for  the  ith  final  user  in  period  t. 

•frit 

= 

the  weight  which  is  applicable  to  the  negative  goal 
discrepancy  for  the  1 final  user  in  period  t. 

0±(t) 

= 

requirement  (goal  amount)  of  support  in  terms  of  its 
dollar  budget  for  the  ith  final  user  in  period  t. 

p^t) 

= 

dollars  budget  of  the  jt^1  producer  in  period  to. 

R 

= 

the  matrix  of  rators  of  support  producers  provides 
to  final  users. 

fT 

= 

vector  of  manpower  per  dollars  of  output. 

The  revised  model  can  be  depicted  as  follows  in  its  transformed  and 
reduced  form  as: 

(Objective  - Minimize  weighted  sum  of  deviations  from  final  user  support 
requirements  and  weighted  sum  of  deviations  from  gross  manpower  goals  and 
weighted  number  of  hires  and  fires  (Reductions-In-Force) ) 

Min  Z Z Q>±t  F+(t)  + *lt  FJ(t))  + Z Z (akt  E^(t)  + 

Bkt  Ek(t)  + ^kt  y(t>  = 6kt  z(t)J 

Subject  to: 

(Final  User  Goal  Constraints) 

Vl(t)  - F+(t)  + F-(t)  = UA(t) 
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(Final  User-Producer  Balancing  Conditions) 
v (t)  = R w (t) 

(Producer  Budget  Constraints) 

Vc> 

(Producer  Manpower  Balancing  Conditions) 
fT  Vj(t)  - tkxk(t) 

(Manpower  Goal  Constraints) 

^(t)  - E^(t)  + E~ (t) 

(Manpower  Transition  Conditions) 


eT  x ( t ) 

(Salary  Budget  Constraints) 
sc  x(t) 

(Non-negativity  Constraints) 


- pj(t) 


- V1’ 


- M x(t-l)  + x(t)  - y(t)  + z(t)  = 0 

(Manpower  Ceiling  Constraints) 


- B(t) 


V1(t),  Wj(t),  xk(t),  yk(t),  zk(t), 

F+(t),  F“(t),  E+(t),  E“(t)  - 0 

This  model  was  extended  further  by  Chames,  Cooper,  Niehaus,  and 
Sholtz  Q- 8.57  to  include  a chance-constrained  form,  which  will  be  discussed 
as  part  of  the  final  section  of  this  paper.  Also  a very  useful  numerical 
prototype  was  developed  £217  which  showed  how  the  models  could  be  integrated 
into  a total  military-civilian  manpower  planning  system.  This  was  a case 
where  the  process  of  developing  the  prototype  itself  helped  to  make  the 
issues  clearer,  and  even  without  further  implementation  efforts  the  results 


had  practical  significance. 
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A series  of  implementation  efforts  were  now  underway.  The  strategy, 
methodology  and  status  of  this  work  was  provided  in  1972  by  Niehaus  . 

Also,  most  of  the  reports  of  the  work  from  its  inception  were  drawn  to- 
gether in  a monograph  by  A.  Charnes,  W.  W.  Cooper  and  R.  J.  Niehaus  /187  and 
also  summarized  in  the  Naval  Research  Logistics  Quarterly  /I? 'J . Implement- 
ation efforts  continued,  with  emphasis  on  local  naval  installations.  These 
included  the  work  at  the  Naval  Underwater  Systems  Center  discussed  by 
A.  Charnes,  W.  W.  Cooper,  R.  J.  Niehaus  and  K.  Padalino  [27j  and  at  the 
Naval  Air  Rework  Facility  at  San  Diego  reported  by  E.  Bres  and  R.  J.  Niehaus 
f4J.  Also  included  was  the  investigation  of  conversational  computer- 
assisted  techniques  by  R.  J.  Niehaus,  D.  Sholtz  and  G.  L.  Thompson  £$lj . 

This  latter  work  grew  into  the  Shore  Activity  Manpower  Planning  System 
(SAMPS)  project  which  is  aimed  at  investigating  large-scale  versions  of 
the  models  at  local  installations  using  the  latest  technology  of  computers 
and  telecommunications  networks » as  discussed  by  R.  J.  Niehaus  and  D. 

Sholtz  Z3l7>  and  R.  J.  Niehaus  /3l7* 

Returning  to  the  mathematical  developments,  parallel  efforts  were  started 
to  find  ways  to  increase  the  computational  efficiency  of  the  solution 
methodology.  One  of  these  efforts  by  A.  Charnes,  W.  W.  Cooper,  D.  Klingman 
and  R.  J.  Niehaus  D-Z7 > D-3J  set  forth  a method  to  calculate  explicit 
solutions  in  convex  goal  programming.  Since  these  references  are  recent 
and  the  ideas  were  also  summarized  in  Z3.&7,  the  comments  in  this  paper  will 
be  limited  to  the  fact  that  the  methods  are  of  much  greater  generality  than 
just  to  the  field  of  manpower  planning.  In  fact,  a non-exhaust ive  list 
of  goal  programs  to  which  these  methodologies  apply  includes: 


a.  Absolute  value  functions,  including  those  with  asymmetric 
weights  and  multi-goal  components; 

b.  General  piecewise  linear  functionals; 

c.  Goal  interval  functions; 

d.  Hypermedium  functionals  and  related  functionals  in  extensions 
of  ordinary  goal  programming. 

The  second  parallel  track  was  to  investigate  methods  using  the  latest 
versions  of  commercial  linear  programming  computer  codes.  The  ideas  from 
the  explicit  solution  research  indicated  that  one  might  be  able  to  secure 
a good  advanced  start  by  the  simple  assumption  that  the  final  solution  will 
satisfy  the  goals.  This  in  fact,  turned  out  to  be  the  case,  and  for  larger 
problems  an  improvement  by  a factor  of  ten  was  generally  possible.  Also, 
it  was  suggested  that  the  use  of  a reduced  dual  form  of  the  model,  which 
has  been  known  for  a long  time  (see  Charnes  and  Cooper  E%J)  , might  secure 
an  added  improvement.—^  This  turned  out  to  be  particularly  true  with  the 
latest  generation  of  commercial  linear  programming  codes  which  use  a 
threaded  list  to  minimize  storage  and  access  time  of  data  during  processing. 

These  developments  obviated  the  need  to  proceed  any  further  with  the 
implementation  of  the  ideas  obtained  in  the  explicit  solution  research.  This 

| 

is  particularly  true,  since  almost  a complete  revision  of  the  input  side  of 
the  computer  support  system  would  have  been  required  to  test  the  explicit 
solution  methodologies  with  an  operational  size  problem. 


2/  See  also  Armstrong  and  Hultz  C"2.J . The  work  on  the  advanced  basis  was 
accomplished  by  S.  Korn  and  the  suggestions  to  use  the  reduced  dual  formu- 
lation were  provided  by  J.J.H.  Forrest  and  J.C.  Jennings  of  Scicon  Computer 
Services,  Ltd.,  during  the  installation  of  the  model  on  the  Computer  Sciences 
Corporation  INFONET  telecommunications  network. 


I 

A 
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In  1975  work  started  on  an  extension  to  the  model  to  incorporate 
equal  employment  opportunity  (EEO)  planning  into  the  goal  programming 
framework.  A comprehensive  model  was  developed  by  A.  Charnes , W.W.  Cooper, 
K.A.  Lewis  and  R.J.  Niehaus  which  included  two  essential  new  features: 

(1)  a set  of  dual  goal  constraints  to  include  both  workload  goals  and  EEO 
goals  by  race-sex  (ethnosexual)  categories;  and  (2)  a "flexibility"  feature 
to  allow  the  model  to  recommend  ways  to  adjust  the  internal  transition 
matrix.  Thus,  both  the  outside  as  well  as  the  inside  possibilities  could 
be  evaluated  through  use  of  the  same  model.  The  definitions  used  in  this 
model  are  as  follows: 

= the  weight  applicable  to  the  personnel  of  ethnosexual 
type  k in  job  category  i in  period  t for  positive 
goal  deviations. 

= the  weight  applicable  to  the  personnel  of  ethnosexual 
type  k in  job  category  i in  period  t for  negative 
goal  deviations. 

Di^(t)  , Dj^t)  = positive  or  negative  deviation  respectively, 
for  the  personnel  of  ethnosexual  type  k in 
job  category  i in  period  t. 


Yikt 

6ikt 

ait 

6it 


the  weight  applicable  to  hires  of  personnel  of  ethno- 
sexual type  k in  job  category  i in  period  t. 

the  weight  applicable  to  fires  (Reductions-In-Force 
(RIF’s))  of  personnel  of  ethnosexual  type  k in  job 
category  i in  period  t. 

the  weight  applicable  to  the  personnel  of  job  category 
in  period  t (associated  with  workload  goals)  for 
positive  deviations. 

the  weight  applicable  to  the  personnel  of  job  category 
i in  period  t (associated  with  workload  als)  for 
negative  deviations. 
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E^(t)  = positive  or  negative  deviations,  respectively, 

for  the  personnel  of  job  category  i (associated 
with  workload  goals)  in  period  t. 

= EEO  goal  for  personnel  of  ethnosexual  type  k in  job 
category  i in  period  t. 

= workload  goal  for  personnel  of  job  category  i in 
period  t. 

= number  of  personnel  of  ethnosexual  type  k in  job 
category  i in  period  t. 

= number  of  hires  of  ethnosexual  type  k in  job  category 
i in  period  t. 

= number  of  fires  (Reductions  in  Force  (RIF's))  of  ethno- 
sexual type  k in  job  category  i in  period  t. 

= initial  inventory  of  personnel  in  ethnosexual  type  k 
in  job  category  i. 

= current  or  "historical"  transition  rate  from  category 
j to  category  i. 

= number  of  personnel  of  type  k in  category  j of  period 
t additionally  transferred  to  category  i for  period  t. 

= number  of  personnel  to  type  k in  category  j not  transiting 
to  category  i in  period  t via  expected  transition  rate 


= "policy  parameter"  to  limit  the  amount  of  additional 
internal  transfers  to  the  "historical"  transition  rate 
in  period  t. 

= minimum  proportion  of  ethnosexual  type  k in  job  category 
i in  period  t. 

= mean  salary  of  job  category  i in  period  t. 

= transfer  costs  (salary  plus  training)  for  the  flexible 
transfers  from  job  category  j to  job  category  i in 
period  t. 

= recruiting  costs  for  new  hires  into  job  category  i in 
period  t. 

= firing  (Reduction  In  Force  (RIF))  costs  from  job  category 
i in  period  t. 
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b*(t) 

total  salary  budget  in  period  t. 

b2(t) 

» 

total  flexible  transfer  (upward  mobility)  budget  in 
period  t. 

b3(t) 

m 

total  hiring  budget  in  period  t. 

b4(t) 

s 

total  f iring(Reduction  In  Force  (RIF))  budget  in 
period  t. 

eT 

s 

sum  vector  (i.e.  (1,1... 1) 

C(t) 

st 

total  manpower  ceiling  in  period  t. 

The  model 

can 

be  written  in  its  transformed  and  reduced  form  as: 

(Objective  - Minimize  weighted  sum  of  deviations  from  EEO  goals  by  ethno- 

sexual  category  including  weighted  hires  and  fires  (Reductions  In  Force) 

and  weighted  sum  of  deviations  from  workload  goals  by  job  category) 

Min  ^ | S £>lkt  D^(t)  + qkt  Dlk(t)  + Yikt  yj(t)  + *J(t)  J 

+ n fa.t  E+(t)  + B ET(t)  J 
it  it  i it  1 

(EEO  Goal  Constraints) 

Xi(t)  ' Dik(t)  + °ik(t)  = Hi(t) 

(Workload  Goal  Constraints) 

l x£(t)  - E+(t)  + E~(t)  = G£(t) 

(Transition  Conditions) 

xi«»  for  all  i,  k = ai 

- j fa±i  Xj(t-1)7  - J q±j  (t)  + 2 r^t) 

+ x^(t)  - y£(t)  + z^(t)  =*  0 

(Maximum  Additive  Flexibility) 

- qjJ(t)  + h^(t)  /J  m117  xj(t-l)  ^ 0 
(Maximum  Subtractive  Flexibility) 

- rjj^t)  + Xj(t-l)  ~ 0 

13 


0 


(Additive-Subtractive  Balance  Conditions) 


J qlj(t)  " i rij(t) 
(Minimum  EEO  Proportions) 

*J(t)  - p*(t)  H*(t) 

Budget  Constraints 

i { 

l ] J IliW 

E E cht)  y*(t) 
i k 1 1 

E E c4(t)  zht) 
i k 1 x 


- 0 

- ^(t) 

- b2(t) 

- b3(t) 

* b4(t) 


(Manpower  Ceiling  Constraints) 

eT  x*(t)  - C(t) 

(Non-negativity  Constraints) 

Dik(t),  D^(t),  E|(t),  E"(t),  x^(t),  q^(t),  r1J(t), 
yk(t)  , z^(t)  are  non-negative  for  all  i,  j,  k,  t. 

The  model  was  then  tested  via  a numerical  example  and  found  to  be 
computable.  Subsequently,  as  outlined  in  Burroughs  and  Niehaus  £ (J , at 
the  request  of  the  Assistant  Secretary  of  the  Navy  (Manpower  and  Reserve 
Affairs) , a reduced  version  of  the  model  was  developed.  The  results  showed 
with  little  doubt  that  the  existing  Navy  civilian  EEO  goals  policy  needed 
substantial  revision.  The  implemented  version  eliminated  the  flexibility 
features  and  also  made  it  possible  to  use  much  of  the  existing  large-scale 
software  system  which  is  in  place.  Rough  estimates  were  made  of  the  EEO 
goals  as  reflected  in  existing  policy.  Essentually,  all  that  was  added  to 
the  Recruiting  Requirements  Model  was  the  EEO  goals  equations,  with  the 
budgetary  constraints  also  eliminated.  Using  the  previous  definitions,  this 
model  can  be  stated  as  follows: 
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(Objective  - Minimize  weighted  sum  of  deviations  from  the  EEO  goals  by 
ethnosexua.l  category  including  weighted  hires  and  fires  (Reductions  In 
Force)  and  weighted  sum  of  workload  goals  by  job  category) 

Mitl  i k t ^ikt  Dlk<t)  + 5ikt  Dik<t)  + Ylkt  yi(t>  + 6ikt  ^(t)7 
+ {I  £*±t  Ei(t)  + Blt  E l(t)J 
(EEO  Goal  Constraints) 

xk(t)  - D.+(t)  + Dt-(t)  - Hk(t) 

(Workload  Goal  Constraints) 

l xk(t)  - E+(t)  + E"(t)  = Gt(t) 

(Transition  Conditions) 

(0)  for  all  i,  k = a^ 

- £ ^mij  xk(t-l)J  + xk(t)  - yk(t)  + zk(t)  = 0 

(Minimum  EEO  Proportions) 

x£(t)  - pk(t)  Hk(t)  - 0 

Manpower  Ceiling  Constraints 

eT  xk(t)  - C(t) 

Non-negativity  Constraints) 

D^t),  Dik(t),  E+(t),  E~ (t)  , xk(t), 

yk(t),  z^(t)  are  non-negative  for  all  i,  j,  k,  t. 

There  was  clear  interest  in  the  Navy  to  move  forward  towards  implement- 
ation. Preliminary  work  was  accomplished  by  Burroughs,  Korn,  Lewis,  and 
Niehaus  Z3  J to  develop  a methodology  to  estimate  the  goals.  The  work  was 
continued  by  Lewis  £28 \J  into  the  development  of  comprehensive  prototypes 
and  of  a systems  concept  for  implementation.  Along  the  way,  it  was  found 
that  the  model  would  not  work  at  the  local  level  due  to  the  small  cell  sizes 
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required.  The  led  to  the  development  of  a coherence  model  by  Charnes, 
Cooper,  Lewis  and  Niehaus  (ViJ  which  uses  "artifact"  goals  in  a capacitated 
distribution  format. 

The  current  implementation  stategy  is  to  develop  a computer  support 
system  for  the  version  without  flexibility  of  the  aggregate  EEO  model  for 
headquarters  planning.  Included  will  be  an  accountability  system  to  put 
the  top  line  officials  of  the  Navy  in  a strong  position  to  monitor  and 
control  progress  towards  the  goals-.  Emphasis  will  also  be  placed  on  the 
development  of  a labor  market  analysis  system  for  use  throughout  the  Navy. 
Research  is  continuing  on  the  models  for  local  installation  planning.  This 
work  will  be  described  in  the  next  section  of  this  report. 

A spinoff  of  the  EEO  prototype  work  was  the  development  of  a promotion 
planning  model.  Here,  the  ethnosexual  categories  were  eliminated  with  a 
model  study  undertaken  to  evaluate  the  promotion  policies  of  the  naval 
laboratory  system  as  described  by  Albanese,  Korn,  Niehaus  and  Padalino  J. 
In  turn,  further  work  is  now  underway  by  Cooper,  Niehaus  and  Nitterhouse 
develop  a better  workforce  goals  planning  system  for  the  naval  laboratories. 
This  is  also  tied  into  the  research  into  conversational  models  in  a tele- 
communications environment  underway  by  Niehaus  Z3l7. 

In  addition  to  the  aggregate  planning  model  applications  selected  for 
feasibility  testing,  there  have  been  a number  of  mathematical  developments 
which  hold  promise  for  the  future.  These  ideas  will  be  discussed  after 
the  next  section  concerned  with  model  of  the  assignment/distribution  type. 


III.  Assignment/Distribution  Model  Structures 


As  was  mentioned  at  the  beginning  of  this  report,  models  of  the 
assignment/distribution  type  have  been  an  important  part  of  the  research 
program.  The  initial  work  by  Charnes,  Cooper,  Niehaus  and  Stedry  ZJ8.£7 
showed  that  the  classical  assignment  model  might  be  extended  in  two  ways. 
The  first  involves  the  use  of  multiple  characteristics  for  each  of  the 
individuals  and  jobs.  The  second  involves  the  dynamic  or  assignment- 
reassignment  aspects  of  career  planning  and  organization  design.  Research 
has  also  been  accomplished  on  ways  mathematically  to  relate  this  planning 
for  individual  persons  and  jobs  with  the  type  of  aggregate  planning  des- 
cribed in  the  previous  section  of  this  report. 

In  order  to  begin  development  of  the  assignment  models,  the  decision 
was  made  to  limit  the  first  numerical  examples  to  the  static  (one-period) 


case.  This  was  done  to  ensure  that  the  solution  algorithms  would  reflect 
properly  the  multiple  characteristic  feature  of  assignment  planning.  At 


this  point,  a biased  quadratic  form  of  the  static  model  was  formulated  by 
Charnes,  Cooper,  Niehaus  and  Sholtz  Z3-8.£7  along  with  a spectral  analysis 
model  for  career  planning.  Following  this  work,  a solution  algorithm  was 
formulated  by  Charnes,  Cooper,  Klingman  and  Niehaus  0\.\] . A statement  of 
this  algorithm  follows: 


part  of  individual  s assigned  to  job  i. 
amount  of  j*"^1  attribute  possessed  by  individual  s. 
amount  of  j1"^1  attribute  desired  in  job  i. 
amount  of  jth  attribute  required  in  job  i. 


weight  on  discrepancy  for  attribute  j in  job  i. 


cis  = "cost"  for  assigning  individual  s to  job  i. 

= set  of  all  attributes  j relevant  to  job  i. 

The  model  can  then  be  stated  as: 

(Objective  - Minimize  the  sum  of  "costs"  across  all  of  the  attributes 

for  all  assignments  made) 

Min  £ E c.  x. 
is  is  1S 

where  c,  = £ a . /a  . - (2r . . + k..)7 

is  jeJ.  SJ  « 

J i 


subject  to 


E x.  =1  Y.  . E x =1  Y (x  = 0,  1 Yi,  s) 


ij  i is 


s is 


0 - r . . - 10  j 0 - agj  - 10  j k±J  - 10  - r±j 


x.  = 0 if  a a . < m. . for  jeJ. 
rs  sj  ij  i 

Two  field  tests  were  conducted  at  naval  installations  as  described  by 
Moore  and  Sholtz  £3Q7  and  by  Bres,  Leader,  Moore  and  Sholtz  C 3.7  • An 
extensive  investigation  of  task  analysis  techniques  was  made  by  Moore  £297 • 
In  the  computer  support  system,  the  very  efficient  capacitated  distribution 
codes  developed  by  Glover,  Karney,  Klingman  and  Napier  £267  were  used. 

The  results  from  the  tests  showed  that  computationally  the  models  could  be 
supported  efficiently.  However,  because  of  the  number  of  attributes 
involved,  in  both  field  tests  the  collection  of  data  was  time  consuming. 
Also,  there  was  a very  apparent  problem  in  trying  to  get  first  line  super- 


visors to  designate  needed  positions:  most  tended  either  to  describe  the 


current  organization,  or  simply  to  furnish  a description  for  each  kind  of 
position  as  formally  defined  in  the  civil  service.  Beyond  these  specific 
problems,  there  was  little  if  any  management  interest  in  the  recommended 
assignments.  In  both  cases  management  could  find  uses  for  by-products  of 
an  assignment  model  based  on  task  inventories,  but  external  constraints  such 
as  unions  and  government  regulations  seem  to  render  less  than  valuable 
the  use  of  assignment  models  for  civil  service  workforces.  This  may  v£ry 
well  not  be  the  case  for  military  workforces,  where  the  assignment  systems 
are  generally  much  more  centralized. 

In  addition  to  the  civilian  personnel  planning,  help  is  being  provided 
to  the  Bureau  of  Naval  Personnel  in  the  design  of  models  for  officer 

distribution.  The  initial  models  were  developed  by  Cass,  Charnes,  Cooper 

/ 

and  Niehaus  L3  3 • Two  models  were  proposed:  (1)  a static  goal  programming 

model  which  is  transformed  into  a distribution  model  for  solution;  and 
(2)  a dynamic  multi-page  model  for  evaluating  officer  rotations.  The  static 
model  has  been  implemented  and  development  is  continuing  on  conversational 
and  bargaining  assignment  extensions  by  Eubanks  and  Thompson  . Further 

assistance  is  being  provided  by  Charnes  and  Cooper  in  continuing  the  develop- 
ment of  the  dynamic  models.  More  recently,  they  are  also  examining  some 
aspects  of  the  enlisted  planning  process. 

During  the  development  of  the  aggregate  EEO  models,  it  became  clear 
that  they  would  not  be  applicable  to  local  installation  planning.  As 
mentioned  previously,  this  was  due  to  the  small  cell  sizes  which  would  be 
required.  However,  from  the  Navy's  management  viewpoint,  it  did  not 
appear  useful  to  develop  upward  mobility  plans  at  headquarters.  It  is  felt 
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that  such  upward  mobility  plans  should  be  developed  locally,  since  that 
is  where  the  civilian  personnel  promotion  and  intra-organizational  transfer 
decisions  are  made.  As  a result,  Charnes,  Cooper,  Lewis  and  Niehaus  Z117 
developed  the  coherence  model,  or  multi-level  EEO  model  (MEEO) . This  model 
is  designed  to  be  in  coherence  with  the  aggregate  headquarters  or  master 
EEO  model.  In  the  MEEO  the  Markov  transition  constraints  are  approximated 
by  "goal  artifacts"  which  relax  the  constraints  to  goals  with  convex  goal 
functionals  on  certain  dyadic  cell  variables.  The  flexibility  features 
which  are  desired  are  also  included  by  means  of  features  of  a form  of  a 
warehousing  model. 

For  this  initial  version  each  ethnosexual  category  is  handled  individually 
by  the  model.  Thus,  to  accommodate  all  the  ethnosexual  categories,  each 
has  to  be  computed  separately  and  added  up  to  obtain  the  organizations' 
manpower  totals.  Later  versions  will  allow  all  the  ethnosexual  categories 
to  be  evaluated  via  a single  model. 

The  model  was  found  to  be  computationally  feasible  on  a small  numerical 
example  in  /1 57  and  in  a more  extensive  example  by  Lewis  /2d7  using  a 
general  purpose  linear  programming  code  for  solution.  With  these  initial 
computational  results  in  hand,  a project  has  been  started  to  build  enough 

of  a computer  support  capability  to  begin  testing  at  a field  site.  The 

preliminary  step  was  to  begin  testing  with  large-scale  highly  efficient 
network  codes  such  as  PNET  £Z6_7,  which  can  handle  multiple  arcs  between 

nodes  as  well  as  lower  and  upper  bounds  on  arc  flows. 

A new,  more  general,  non-linear  goal-arc  network  model  has  been  developed 
by  Charnes,  Cooper,  Nelson  and  Niehaus  /j.7 'J  to  be  used  in  connection  with 
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the  PNET  code.  This  model  which  is  an  extension  of  the  version  suggested 
by  Charnes,  Cooper,  Lewis  and  Niehaus  uses  the  following  variables: 


x..(t)  = the  number  of  personnel  of  type  a transferred 

from  job  category  i to  job  category  j in  period  t; 

= the  current  or  "historical  transition  rate"  from 
job  category  i to  job  category  j; 

p^  = the  desired  or  actual  proportion  of  personnel  of 
type  a in  job  category  i; 

a^t)  = the  aggregate  workforce  goal  for  job  category  i in 
period  t; 

g^j(t-l,t)  = a^(t-l)  ^ = the  number  of  personnel 

of  type  a expected  to  transfer  from  job  category  i 
to  job  category  j in  period  t; 

p^a^(t)  = the  projected  workforce  goal  for  personnel  of  type 

a in  job  category  i in  period  t; 


xoi(l)  = the  number  of  personnel  of  type  a initially 

aboard  in  job  category  i; 

iQ  = subscript  indicating  natural  attrition 

jQ  = subscript  indicating  outside  source  or  firing; 

y^(t)  = the  total  number  of  personnel  of  type  « in  job 

category  i in  period  t. 

To  each  job  category  in  each  period  we  assign  two  nodes  an  antecedent 
and  a consequent.  We  designate  the  class  of  antecedent  "job"  nodes  for  period 
t as  J~(t);  the  class  of  consequent  "job"  nodes  by  .^(t).  Jj“(t)  is  the  ith 
job  antecedent  node;  Jj+(t)  is  the  j11*1  job  consequent  node. 

Xjj(t)  = the  flow  from  nodeJj~(t)  to  node  Ji+(t)  on  the  k^  individual 
arc  of  the  "goal  arc". 
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For  each  proper  (real)  job  between  two  periods  we  designate  a 
"value"  node  to  receive  the  goal  arc  flow  from  the  consequent  node  of  the 
immediate  past  period.  Let  V^(t)  denote  the  "value"  node  for  job  i between 
period  t-1  and  t. 

yk(t)  = the  flow  on  arc  k from  node  J^+(t-l)to  node  V^(t). 

L^(t)  and  U^t)  are  lower  and  upper  on  the  flow  from  Jj“(t) 

to  Jj+(t)  via  the  kth  arc.  For  some  k,  U^(t)  = g^t-l.t).  L^(t)  and 

k 

U^t)  are  lower  and  upper  bounds  respectively  on  the  flow  from  J^+(t)  to 
the  "value"  node  V^t).  From  some  k,  U^(t)  = p^a^(t)  . 

The  model  can  be  written  as: 


(Objective  - minimize  the  sum  of  deviations  of  actual  transitions  from 
expected  transitions  between  job  categories  over  all  times  periods  and 
the  sum  of  deviations  of  the  number  of  personnel  in  job  categories  from 
the  targeted  workforce  goal  for  the  job  categories  over  all  time  periods.) 


Min  Z E Cj  x . ( t)  + E £ E 


t k j iH0  _1  “iJ 


k t i^ic 


4 y£(t) 


Subject  to 
(Goal  Conditions) 


1 x^t)  - xJ4(t)  = 0 


ij 


l ^i(t)  - ^i(t) 


= 0 


(Manpower  Flow  Conditions) 


E,  £ Xji(n)  + E x..  (t)  - E x . (1)  = 0 

j^iQ  i ij  t JV  ieJ-(l)  oi 


soi 


jtJ+(l)  k 


E xjd) 


E E 

k iej-(l) 


x^d)  - E yf(t)  = 0 


where  j ^ j 
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One  of  the  next  steps  will  be  to  bring  together  the  aggregate  planning 
and  assignment/distribution  type  models  into  one  model  system.  Such  systems 
have  been  suggested  by  Cass,  Charnes , Cooper  and  Niehaus  C ~J  for  officer 
distribution  planning  and  by  Charnes,  Cooper,  Lewis  and  Niehaus  ZI&7  for 
EEO  planning. 


IV.  Other  Extensions  and  Research  Frontiers 

In  a long  term  modeling  research  program  such  as  this  one,  many  other 
extensions  have  been  developed  without  implementation,  full  or  partial. 

This  has  been  due  to  either  lack  of  resources  for  implementation,  change 
in  focus  of  the  project,  or  the  need  for  more  mathematical  development  than 
could  be  accomplished  at  the  time.  Several  of  these  ideas  will  be  reviewed 
so  that  a discussion  of  the  research  frontiers  might  be  more  complete. 
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The  idea  of  planning  under  uncertainty  or  risk  was  advanced  in  1970 
by  Charnes,  Cooper  and  Niehaus  /l8.^7  and  developed  further  by  Charnes, 
Cooper,  Niehaus  and  Sholtz  /18.57  through  the  introduction  of  the  use  of 
chance-constrained  programming.  In  particular,  a discussion  was  provided 
in  £8.#  of  extending  the  Markov  matrix  notations  of  the  Recruiting 
Requirements  Model  to  chance-contrained  intrepretations . In  £18.^7  these 
chance  constrained  applications  were  extended  further  to  deal  with  risk 
variations  in  the  right  hand  sides. 

One  means  of  doing  this  is  by  employing  "zero-order  decision  rules"  of 
chance-constrained  programming.  Under  this  rule  the  values  of  the  decision 
variables  are  all  chosen  in  advance  of  knowledge  of  the  sample  values  of 
the  random  variables.  Using  the  notation  of  the  Recruiting  Requirements 
Model,  we  can  examine  the  new  elements  which  the  use  of  chance-constrained 
programming  introduces.  For  illustration,  the  k11*1  Manpower  Goal  would  be 
written  as: 

p rxk(t)  ± Gk(t)7  * a)kt 

This  means  that  the  "kt^"  manpower  on-board  in  period  t must  be  chosen  so 
that  it  does  not  exceed  the  kC^  Manpower  Goal  in  period  t with  a probability 
of  at  least  «kt.  Note  that  Gk(t)  is  a random  variable  whose  sample  value  is 
not  known  when  the  planning  decision  for  xk(t)  must  be  made.  It  is  required, 
however,  that  x^(t),  when  selected,  must  not  exceed  Gk(t)  by  a probability 
<okt  which  is  also  stipulated  prior  to  the  knowledge  of  the  sample  value 
of  G^t).  Only  the  probability  distribution  for  Gk(t)  is  known  when  the 
planning  decisions  are  to  be  made. 
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With  the  zero-order  decision  rules,  the  chance  constraints  may  be 
inverted  to  obtain  a new  linear  programming  problem  of  virtually  the  same 
structure  as  the  deterministic  case.  The  details  of  this  process  can  be 
found  in  Zl8.57.  The  testing  of  the  chance-constrained  extension  was  not 
accomplished  because  of  more  immediate  implementation  considerations.  Now, 
however,  the  civilian  manpower  modeling  program  is  again  moving  toward  the 
frontiers  of  goal  setting  in  the  areas  of  EEO  and  program  planning.  At 
an  appropriate  point,  at  least  a numerical  prototype  will  be  constructed 
to  check  such  issues  as  management  acceptance,  data  availability,  and 
computer  support  requirements. 

Tracing  another  thread,  it  is  interesting  to  see  how  an  idea,  dormant 
for  several  years,  reappears  in  an  improved  form.  Such  is  the  case  with 
the  possibility  observed  in  1971  by  Charnes,  Cooper,  and  Niehaus  /!207  of 
redesigning  the  system  transition  matrices  to  obtain  a stable  desired  mix 
of  personnel.  In  this  case  one  wishes  to  insure  "as  closely  as  possible" 
the  possibility  of  effecting  transfer  and  the  costs  of  influencing  transfer 
rates  subject  to  constraints  on  available  resources  to  obtain  this  stable 
mix.  This  was  rephrased  mathematically  by  letting  the  vector  of  desired 
proportions  be  represented  by  n and  the  desired  matrix  of  transitions  by  1 
Ideally,  for  a steady  state,  or  equilibrium  result,  this  would  be 
7T  = Mn 

If,  in  any  actual  situation,  the  current  transition  matrix  is  designated 
by  M0,  the  possible  changes  will  be 
M = M - M 

o 
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with  the  constraint  system 
G (AM)  - H 

where  G is  a coefficient  matrix  and  H is  a matrix  of  constants. 

The  above  equations  then  lead  to  the  goal  programming  problem: 

(Objective  - minimize  the  weighted  discrepancies  between  the  desired 
proportions  and  the  desired  proportions  multiplied  by  the  system  transition 
matrix) 

Min  w^  | tt  - Mn  | 

subject  to 
(Constraint  System) 

G (AM)  ^ H 

This  model  can  be  seen,  of  course,  as  the  initial  conceptual  development 
of  the  flexibility  features  which  are  included  in  the  EEO  models  and  in 
the  promotion  planning  model. 

The  above  example  illustrates  the  desirability  of  maintaining  a funda- 
mental mathematical  research  program  as  part  of  an  on-going  system  develop- 
ment which  has  many  parts  operational  or  nearing  implementation.  There  are 
still  many  other  frontiers  which  are  worth  exploring  for  future  improvements. 
The  work  on  explicit  solutions  to  convex  goal  programs  will  more  than  likely 
reappear  as  larger  and  larger  systems  are  brought  together.  Also  on  the 
horizon  might  be  the  integration  of  information  theoretic  statistics  with 
mathematical  programming  as  suggested  by  Charnes  and  Cooper  [})]  in  their 
other  research. 

Still  other  scientific  disciplines  are  being  integrated  into  the  Navy 
Civilian  manpower  modeling  research  program.  It  is  clear  that  the  fields 


of  behavioral  science,  managerial  economics,  management  information,  and 
labor  economics  are  as  much  a part  of  the  effort  as  is  the  mathematical 
research  on  which  this  paper  has  concentrated.  These  interactions  might 
better  be  held  as  subjects  of  continuing  reports.  It  is  fitting,  however, 
to  conclude  this  paper  by  saying  that  the  pioneering  efforts  which  make 
this  integration  possible  would  not  have  been  accomplished  without  the 
contributions  of  Abe  Charnes  and  Bill  Cooper.  In  all  likelihood  their 
efforts  will  continue  to  set  the  stage  for  the  1990's  and  beyond. 
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